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Black holes have been found over a wide range of masses, from stellar remnants with 
masses of 5-150 solar masses (M.,), to those found at the centres of galaxies with 
M>10°M.. However, only a few debated candidate black holes exist between 150M, 
and 10°M.,. Determining the population of these intermediate-mass black holes is an 
important step towards understanding supermassive black hole formation in the 
early universe’”. Several studies have claimed the detection of a central black hole 

in w Centauri, the most massive globular cluster of the Milky Way* >. However, these 
studies have been questioned because of the possible mass contribution of stellar 
mass black holes, their sensitivity to the cluster centre and the lack of fast-moving 
stars above the escape velocity® *. Here we report the observations of seven fast- 
moving stars in the central 3 arcsec (0.08 pc) of w Centauri. The velocities of the fast- 
moving stars are significantly higher than the expected central escape velocity of the 
star cluster, so their presence can be explained only by being bound toa massive black 
hole. From the velocities alone, we can infer a firm lower limit of the black hole mass 
of about 8,200M.,, making this a good case for an intermediate-mass black hole in the 


local universe. 


w Centauri (w Cen) is a special case among the globular clusters of 
the Milky Way. Owing to its high mass, complex stellar populations 
and kinematics, œw Cen is widely accepted as the stripped nucleus of 
an accreted dwarf galaxy’°". These factors combined with its prox- 
imity (D = 5.43 kiloparsec (kpc); ref. 12) have made it a prime target 
for searching for an intermediate-mass black hole (IMBH). As part 
of the oMEGACat project”, we recently constructed an updated 
proper-motion catalogue of the inner regions of w Cen, on the basis of 
more than 500 Hubble Space Telescope (HST) archival images taken 
over atime span of 20 years. The unprecedented depth and precision 
of this catalogue have enabled us to discover a statistically significant 
overdensity of fast-moving stars in the centre of the cluster (Fig. land 
Extended Data Fig. 1). In total, we find seven stars with a total proper 
motion higher than 2.41 mas yr‘ within 3” of the centre determined in 
refs. 6,15 (hereafter AvdM10 centre). At a cluster distance of 5.43 kpc 
(ref. 12), this corresponds to projected two-dimensional (2D) veloci- 
ties higher than the escape velocity of the cluster ifno IMBHis present 
(Vese = 62 km s~”; ref. 16; Methods). Here we show that the presence of 
these stars strongly indicates a massive black hole, similar to the S-stars 
in the Galactic centre”. A list of the fast-moving stars is provided in 
Extended Data Table 1, and we label the fast-moving stars with letters 
from A to G, sorted by their proximity to the AvdM10 centre. All these 


stars lie along the cluster main sequence in the colour-magnitude dia- 
gram (CMD) (Fig. 2). The fastest and the centremost star (Fig. 1, star A) 
has a 2D proper motion of 4.41 + 0.08 mas yr” (113.0 + 1.1 km s”). The 
motion of this star was measured over 286 epochs and a full 20.6-year 
time baseline (Fig. 3). We run extensive quality checks to ensure that the 
astrometry of the discovered fast-moving stars is reliable. To ensure the 
cleanest possible dataset, we limit our analysis to stars whose velocity 
is at least 3g above the escape velocity. This leads to the exclusion of 
stars Band G; however, this has negligible influence on the determined 
IMBH constraints. 

Four of the fast-moving stars, including the three fastest in the 
sample, are found within the centremost arcsecond (Fprojectea < 0.03 pc 
OF Fprojected < 0-09 ly). These four innermost stars are all fainter than 
Meeogw > 22.7, which is unlikely (P= 0.013) to be a random occurrence 
given the overall distribution of stellar magnitudes in the centre of 
w Cen. Moreover, all of them lie towards the blue side of the main 
sequence. Both these properties could have interesting physical impli- 
cations on the mechanism involved in capturing these stars or on their 
tidal interactions with the IMBH. 

We expect a certain number of Milky Way stars in our field of view, 
and because they have a large proper motion with respect to w Cen, 
they can mimic fast-moving cluster stars. On the basis of the number 
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Fig. 1|Location of fast-moving stars. a, Solid grey markers indicate all stars 
witha high-quality (HQ) proper motion measurement in our new catalogue 
within a 40” x 40” region centred on the AvdM10 centre®. Well-measured stars 
with velocities higher than the cluster escape velocity (62 km s”) and lying on 
the cluster main sequence in the colour-magnitude diagram (Fig. 2) are marked 
in pink. We use filled markers for stars that are at least 30 above the escape 
velocity, and open markers for all other stars above the escape velocity. 

b, Stacked image of the innermost region of w Cen using all observations in 

the WFC3/UVIS F606W filter. The fast-moving stars and their proper-motion 


density of fast-moving stars at larger radii (Extended Data Fig. 1), we 
estimate the rate of contaminants to be 0.0026 arcsec”, which is con- 
sistent with the expectations from the Besancon Milky Way model”. 
This number density gives an expected average value of only 0.074 
foreground stars in the inner 3-arcsec radius. A detection of five such 
stars by a pure coincidence can, therefore, strongly be ruled out by 
simple Poisson statistics (P = 1.7 x 10°°) (ref. 19). Having two or more 
random contaminants within our five-star sample can also be ruled 
out at the 3g level (P= 0.0026). We also show in the Methods that these 
stars cannot be explained by objects bound to stellar mass (S100M.,) 
black holes and that ejections from three-body interactions or anIMBH 
are not plausible. 

Therefore, the presence of the seven central stars moving faster 
than the escape velocity of the cluster can be explained only if they 
are bound to a compact massive object near the centre, raising the 
local escape velocity. If no massive object was present, their velocities 
would cause them to leave the central region in less than 1,000 years, 
and then eventually escape the cluster. These fast-moving stars area 
predicted consequence of an IMBH, but these stars are not expected 
from mass-segregated stellar-mass black holes’. 

We do not know several parameters of the system, including the 
mass and exact location of this massive object, the relative line-of-sight 
(LOS) distance between the object and the stars, and the LOS velocity 
of the stars. Despite this, we can calculate a lower limit on the mass of 
the dark object using the 2D velocities of the fast-moving stars, assum- 
ing only that they are bound to it. The combined constraint from the 
velocities of the five robustly measured fast-moving stars is about 
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vectors are shown in pink. The arrows indicate the stellar motion over 100 years. 
We also list the measured individual velocities. The cyan cross indicates 

the photometric centre of w Cen measured by ref. 6, with the dashed circle 
indicating the 1” error reported for this centre; the orange cross marks the 
centre, allowing for the lowest IMBH mass; and the blue cross marks the most 
likely position of an IMBH given the Markov chain Monte Carlo (MCMC) analysis 
of the acceleration limits of the fast-moving stars. Dec, declination; RA, right 
ascension, min., minimum. 
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Fig. 2 | HST-based colour-magnitude diagram (CMD) of w Cen. The CMD 
locations of all fast-moving stars are marked with pink symbols, with their 
photometric 1o errors marked with error bars. All of them lie on the main 
sequence, showing that they are probably members of w Cen. The stars are 
labelled from A to G, sorted by their distance from the AvdM10 centre. 
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Fig. 3 | Motion of the fastest star. a, Individual measured positions with 
loerror bars. b, Multi-epoch HST imaging for star A, the fastest (V projected = 
113.0 + 1.1km s”) and centremost of the seven fast-moving stars discovered in 
the centre ofw Cen. The star is indicated witha pink marker on the images and 


8,200M., thus making an IMBH the only plausible solution. The posi- 
tion of the IMBH requiring the lowest mass is only 0.3 arcsec away 
fromthe AvdM10 (ref. 6) centre, in agreement with the +1” error onthe 
AvdM10 centre. Further details of this calculation can be found in the 
Methods. 

Although the linear motion and the velocity of the stars can be meas- 
ured with great precision, the expected acceleration signal from an 
IMBH is considerably weaker and harder to detect. However, even a 
non-detection of acceleration could provide useful constraints on 
the mass and location of the IMBH. The accelerations of all stars are 
consistent with zero within 30, but two stars have more than 2o accelera- 
tion measurements. We model both the velocity and the acceleration 
measurements to further constrain the IMBH properties (Methods); 
this calculation increases the lower limit on the black hole mass to 
21,100 M. (99% confidence) and gives a preferred position for the IMBH 
0.77” northeast of the AvVdM10 centre. 

Apart from these constraints that are purely based on the assumed 
escape Velocity and our astrometric data of the five robustly meas- 
ured fast-moving stars, we also compared the full velocity distribution 
observed in the inner 10” of w Cen to already existing state-of-the-art 
N-body models? with various IMBH masses. Models with no IMBH, 
a stellar mass black hole cluster or an IMBH with a mass greater than 
50,000M, are all strongly ruled out, whereas models with an IMBH 
mass of 39,000M., and 47,000M. are most consistent with the fraction 
of fast-moving stars and the observed velocity distribution. However, 
we caution that our comparisons show that low number statistics 
limit these comparisons, and mismatches with the overall velocity 
distribution suggest a need for improved modelling (see Methods 
for more details). 

The detection of fast-moving stars in the centre of w Cen strength- 
ens the evidence for an IMBH in this cluster. Owing to the probable 
origin of w Cen as the nucleus of the Gaia~Enceladus-Sausage dwarf 
galaxy”°”!, this black hole provides an important data point in the 
study of black hole demographics in low-mass galaxies, along with 
other black holes that have been detected in more massive globular 
clusters and stripped nuclei around M31 such as G1 (M = 20,000M.,) 
(refs. 22,23) or BO23-G078 (M = 100,000M.) (ref. 24). Moreover, this 
black hole provides the closest massive black hole and only the sec- 
ond after Sgr A* for which we can study the motion of multiple indi- 
vidual bound stellar companions. A comparison with the motion of 
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its motion over 21 years with a line. This plot shows the notable astrometric 
quality and the long temporal baseline we have in our unique dataset. Similar 
plots for all other stars are shown in Extended Data Fig. 2. Dec, declination; 
RA, right ascension; FOV, field of view. 


the stars in the Galactic centre is shown in the Methods and Extended 
Data Fig. 8. 

A more precise estimate of the black hole mass requires dynami- 
cal modelling of all newly available kinematic data using models that 
include the impact of both an IMBH and mass-segregated dark rem- 
nants. The exact properties of the orbits of the fast-moving stars have 
to be determined by deep, pin-pointed follow-up observational studies. 
Spectroscopic observations with integral-field-unit instruments such 
as VLT MUSE” or JWST NIRSpec IFU” could yield LOS velocities for the 
fast-moving stars. Even more precise and deeper astrometric measure- 
ments with existing (VLTI GRAVITY+’, JWST NIRCam7’®) or future (ELT 
MICADO”’, VLT MAVIS”°) instruments could enable the detection of 
additional tightly bound stars, and the measurements of accelerations, 
key for obtaining direct measurements of the black hole mass. Our 
results also motivate revisiting the other likely accreted nuclear star 
clusters of the Milky Way”, with M54 being the clearest case. For the 
search for IMBHs in other globular clusters, our results indicate that 
it may be necessary to extend kinematic studies to the faintest stars, 
which is observationally challenging for clusters at larger distances 
and with high central densities. 
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Methods 


Discussion of previous IMBH detections in w Centauri 

The debate about an IMBH in w Cen dates back almost two decades 
but has remained controversial. Early dynamical modelling based 
on LOS integrated-light velocity dispersion measurements suggested 
anIMBH mass of Mypy = (4.0'02°) x 104M, (ref.3). These results were 
challenged with a precise redetermination of the centre of the 
cluster® and dynamical modelling of proper motions’ measured from 
multi-epoch HST imaging observations that placed an upper limit of 
1.2 x 104M, on the IMBH. Using additional integrated light observa- 
tions and a centre based on the maximum LOS velocity dispersion, 
a best-fit IMBH mass of Myypy = (4.741.0) x 10*M, was obtained in 
ref. 4. When assuming the AvdM10 centre, the IMBH mass was slightly 
lower, Mimpu = (3-0 £0.4) x 104M.. 

Subsequent comparisons of both proper motions and LOS velocities 
to N-body simulations continued to show evidence for an approxi- 
mately 4.0 x 104M, IMBH*!2. However, these observations were also 
shown to be fully consistent with a dark cluster of stellar mass black 
holes in the central region of w Cen®. The lack of fast-moving stars in 
previous proper-motion catalogues supported this scenario over an 
IMBH”. Other works noted the influence of radial velocity anisotropy 
on dynamical mass estimates. 

Most recently, the discovery of a counter-rotating core using VLT 
MUSE LOS velocity measurements of individual stars” highlighted 
once again the kinematic complexity of the centremost region of w Cen. 
The centre of this counter-rotation coincides with the AvdM10 centre 
within about 5” but is incompatible with the centres used in refs. 3,4. 


Previous accretion constraints in context 

With the detection of a10*°M, IMBH, the upper limits on any accretion 
signal at X-ray% and radio” wavelengths make this the most weakly 
accreting black hole known. Deep, 291 ksec Chandra observations place 
an upper limit of the 0.5-7 keV luminosity of around 10” ergs s 
(ref. 36), roughly 12 orders of magnitude below the Eddington limit. 
The radio upper limit implies an even fainter source, with the 5 GHz 
upper limit of 1.3 x 10” ergs s” (ref. 37) corresponding to an implied 
X-ray luminosity using the fundamental plane of about 10” ergs s7 
(ref. 38). Assuming standard bolometric corrections of about 10 
(ref. 39), this X-ray luminosity upper limit suggests an Eddington ratio 
of log(Lpo\/Leaa) < - 12, far fainter than that for Sgr A**° or any other 
known black hole. This faint signal could be because of a combination 
of low surrounding gas density, alow accretion rate of that gas and/or 
alow radiative efficiency”. Low-luminosity active galactic nuclei includ- 
ing Sgr A* are brightest at infra-red and sub-millimetre wavelengths 
most likely because of synchrotron emission from compact jets”. 
Therefore, future observations with the James Webb Space Telescope 
or the Atacama Large Millimeter array would provide the highest sen- 
sitivity to any emission from the IMBH of w Cen. Any detection would 
reveal the location of the IMBH as well as provide valuable constraints 
on the black hole accretion in this extremely faint source. 


Proper-motion measurements and sample selection 

Our proper-motion measurements are based on the reduction of archi- 
val HST data of the central region of w Cen, taken over a time span of 
more than 20 years. We used the state-of-the-art photometry tool KS2 
(ref. 32) for the source detection and the astro-photometric measure- 
ments, and the established procedure described in refs. 43-46 to meas- 
ure proper motions relative to the bulk motion of the cluster. The result 
of this extensive study is a proper-motion catalogue with high-precision 
measurements for 1.4 million stars out to the half-light radius of w Cen 
with a typical temporal baseline of more than 20 years. Owing to the 
large number of observations (in total we reduced over 500 images 
and some stars in the central region have up to 467 individual astro- 
metric measurements) the catalogue reaches unprecedented depth 


and precision. The highest precision is achieved in the well-covered 
centre of the cluster, in which our proper motions have a median error 
of only about 6.6 u as yr” (0.17 km s”) per component for bright stars. 

The catalogue is larger than any other kinematic catalogue pub- 
lished for a globular cluster and significantly extends previous proper- 
motion catalogues for w Cen®**”, A detailed comparison with other 
proper-motion datasets is published along with the catalogue”. Ina 
following section and in Extended Data Fig. 3, we compare the com- 
pleteness of the different catalogues to show that it is plausible that 
the fast-moving stars have been missed in previous searches. 

We use a high-quality subset of the proper-motion catalogue to 
search for real fast-moving stars and limit spurious astrometric meas- 
urements (for example, two sources that are falsely identified as one) 
that can have apparent high proper-motion measurements. Our criteria 
for this subset are based on the amount of available data for the meas- 
urements. Specifically, we used only sources that had at least 20 astro- 
metric measurements covering a temporal baseline of at least 20 years 
anda fraction of rejected measurements (based on sigma clipping) of 
less than 15%. We also made cuts on the quality of the proper-motion fit 
requiring both a proper-motion error less than 0.194 mas yr“ = 5 km s” 
andareduced x’ < 10 for the linear proper-motion fit for both the right 
ascension (RA) and declination (Dec) measurements. Apart from these 
quality selections, we also required the star to lie on the CMD sequence 
in an HST-based CMD (Fig. 2). These cuts help to drastically reduce the 
number of contaminants. These criteria are met uniformly out toa 
radius of about 90 arcsec; at larger radii they lead to selection effects 
due to reduced observational coverage. A total of 157,320 out of 241,133 
(65.2%) entries of the proper-motion catalogue within r< 90” match 
the combined criteria. Extended Data Table 1b shows the individual 
measured proper-motion components for the seven fast-moving stars. 
We note that for this analysis we have not applied the local a posteriori 
proper-motion corrections provided with the catalogue”, as we are 
studying the central region that is well dithered and observed with vari- 
ous rotation angles. We verified that applying these corrections would 
neither change our fast-moving star sample nor alter our conclusions. 


Details on verification for fast-moving stars 

The criteria detailed above should lead to a clean dataset with very 
few spurious proper-motion measurements. To ensure that the meas- 
urements for the fast-moving stars are reliable, we inspected each of 
them carefully. 

As afirst step, we tested the quality of the raw astrometric measure- 
ments by studying several goodness-of-fit parameters and photometric 
quality indicators for the point-spread-function fits used to measure 
stellar positions (Extended Data Fig. 3). We performed this analysis for 
the WFC3/UVIS F606W filter as it is the most used filter in the centre of 
w Cen, and each star has at least 195 measurements in this filter. To 
verify the goodness of fit, we used the mean of the so-called quality-of- 
fit (QFIT) flag and the radial excess value, both of which take into 
account the residuals of the point-spread-function fit. Furthermore, 
we looked at the mean of the ratio of source flux with respect to the flux 
of neighbouring sources within its fit aperture. All five stars used for our 
analysis behave typically for well-measured stars of their magnitude and 
none of them show extreme values that would indicate problems with 
the photometry. It is noteworthy that the two stars excluded from our 
analysis based on their velocities being less than 30 above the escape 
velocity show some deviations. Star B has a relatively low mean QFIT 
value and high radial excess. Star G is the only one in the sample for 
which the neighbour flux to star flux ratio is larger than 1. 

Asasecond step, we looked at the stars in several stacked HST images 
taken at various epochs ranging from 2002 to 2023. The extensive 
multi-epoch imaging is demonstrated in Extended Data Fig. 2. The 
proper motion ofa star at the escape velocity (62 km s”) is 2.41 mas yr“. 
Therefore, we expect to see a displacement of at least 50 mas over 
21 years, which corresponds to 1.25 WFC3/UVIS pixels. This motion can 
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be seen by eye for all seven stars in the multi-epoch images. Again the 
excluded stars B and Gstick out, that is, that they are partially blended 
with the neighbouring stars, thus explaining their larger astrometric 
errors. 

Finally, we tested the reliability of our proper-motion measurements 
by limiting the raw position measurements to different subsets and 
redoing both the linear and quadratic fits to the motion of the stars. 
The first test run included only high signal-to-noise measurements. The 
second test runincluded only measurements taken with the WFC3/UVIS 
F606W filter. By using only one filter, we are immune to colour-induced 
effects such as a partially resolved blend between two differently col- 
oured stars. The proper motions of all fast-moving stars are consistent 
within the measurement uncertainties using both methods. 


Comparison with other proper-motion datasets 

Before our analysis, two other high-precision proper-motion cata- 
logues based on HST data have been published®” covering the centre 
of w Cen. Both datasets were searched for central high proper-motion 
stars, but none of the stars in our sample have been reported before. To 
understand why this is the case, we compare the completeness of the 
different catalogues. Extended Data Fig. 3a shows histograms of the 
magnitudes of stars with measured proper motions. Inthe inner 20”, our 
new catalogue contains more than three times the number of stars of the 
literature catalogues and extends to significantly fainter magnitudes. 
The newly detected fast-moving stars all lie at faint magnitudes, in which 
the completeness of the older catalogues is significantly lower thanin 
the new proper-motion catalogue. This is because of the larger amount 
of data and the updated source-finding algorithms in the new catalogue, 
which explains the previous non-detection of the fast-moving stars. 


Discussion of photometric errors 

Apart from the astrometric reliability, we also studied the quality of 
the photometric measurements used to locate the stars in the CMD. 
Although the innermost stars are faint, their statistical photometric 
errors are small because of the large number of individual photometric 
measurements combined to a weighted mean value. The statistical 
errors range from 0.004 mag to 0.037 mag and are given in Extended 
Data Table 1c. However, especially for faint stars, this statistical error 
is not able to capture systematic issues caused, for example, by the 
influence of brighter neighbouring stars. These issues can be identi- 
fied only by verifying the quality of the point-spread-function fit used 
to determine the individual photometric measurements. We report 
the mean quality of fit, radial excess and neighbour flux to source flux 
ratio flags for both filters in Extended Data Table 1c and compare them 
with those of stars at similar (Am < 0.5) magnitudes in Extended Data 
Fig. 3. All stars in the robustly measured sample show typical quality 
of fit for their respective magnitude. We note that stars B and G (which 
were excluded from the analysis) show comparatively poor QFIT. Stars 
E and G show a possible flux contribution from a neighbouring source 
(indicated by a high radial excess value and a high neighbour flux to 
source flux ratio). This can be confirmed by the stacked images shown 
in Extended Data Fig. 2, in which these stars show a close neighbour. 
Owing to the relatively bright magnitude of star E and the low astro- 
metric scatter, we still consider its measurement valid. 


Density of Milky Way contaminants versus fast-star background 
To quantify our expected level of contamination from Milky Way fore- 
ground and background stars, we compared our results with those of 
a Besancon model”. Using the m1612 model, we simulate a 1 square 
degree patch centred on w Cen retrieving Johnson colours and kin- 
ematics. We then transform the model Johnson Vand / magnitudes 
into F606W and F814W magnitudes using linear relations fitted to 
Padova models* between V-/of 0 and 2. We use the same colour cuts 
and consider stars between F606W of 16-24. We then count the number 
of stars with a total proper motion above 2.41 mas yr“ (equivalent to 


our velocity cutoff at the escape velocity Ves = 62 kms‘). These would 
appear as contaminants in our fast-moving star sample. We find a 
density of 0.0039 stars per arcsec”. This is somewhat higher than the 
0.0026 + 0.0003 stars per arcsec? found as the background level in our 
observations. This discrepancy is alleviated by considering only 65.2% 
of allstars within our catalogue meet the high-quality criteria used for 
our fast-moving star selection (see above). Correcting for this factor, 
we get an expected background of 0.0025 stars per arcsec’, perfectly 
matching the observed background density (Extended Data Fig. 1). This 
suggests that our background level is consistent with being predomi- 
nantly Milky Way contaminants. Relaxing our requirement that stars 
be more than 3g above the escape velocity results in a higher observed 
background level of 0.0042 stars per arcsec’, no longer consistent with 
the Milky Way background. This suggests that our stricter definition 
of a fast-moving star reduces contamination from poorly measured 
stars in w Cen to a negligible level. 


Other scenarios that could explain the fast-moving stars 

A complete contamination of our sample by Milky Way foreground 
and background stars that are non-members of w Cen can be ruled out 
statistically. We now explore and rule out alternative scenarios to the 
fast-moving stars being bound to an IMBH. One alternative explanation 
for stars with a high velocity is to have them bound in a close orbit with 
a stellar-mass black hole. This scenario can be ruled out for BHs less 
than 100M,, as the periods required to reach the observed velocities 
are less than 10 years, which is well within the 20-year span over which 
we have observed linear motions. 

Another scenario could be that the stars are actually unbound from 
the cluster and have recently been accelerated by three-body interac- 
tions, either with stellar-mass black hole binaries or with an IMBH. 
Ejection by an IMBH in the centre of the cluster can be ruled out by the 
very high rate of ejections necessary to sustain the observed number 
of fast-moving stars within the centre and the absence of observed 
fast-moving stars at larger radii. To sustain a density of 0.18 fast-moving 
stars per arcsec’ in the inner 3 arcsec (equivalent to our conservative 
sample of five stars) moving with at least 2.4 mas yr“ would require 
ejections with a rate of 0.004 stars per year. This would lead to about 
117 additional fast-moving stars at larger radii (20” < r < 90”), apart 
from around 60 foreground stars expected from the (completeness 
corrected) Besancon Milky Way model. In our dataset, we find 61 
fast-moving stars between 20” < r < 90”, consistent with the expected 
Milky Way background but not consistent with a substantial number 
of additional ejected stars. Moreover, a high hypothetical ejection rate 
of 0.004 stars per year would deplete all of the about 10 million stars of 
w Cen injust 2.5 Gyr. If no IMBH is present, accelerations of stars above 
the escape velocity are still possible by three- or four-body interac- 
tions between stellar or compact object binaries*’. However, these 
interactions would not be limited to the innermost few arcseconds of 
the cluster because of the slowly varying stellar density in the core of 
w Cen. Furthermore, the expected rate of these ejection events is of the 
order of less than one ejection per 1 million years, around 1,000 times 
lower than needed to explain the observed number of fast-moving 
stars in the centre of w Cen*”°. 


Search for the fast-moving stars in recent LOS velocity data 

Line-of-sight (LOS) velocities of the fast-moving stars could help to 
exclude contaminants and provide further constraints on the orbits 
of the stars and the mass and position of the IMBH. The deepest and 
most extensive spectroscopic catalogue of stars in w Cen is part I of 
our recently published oMEGACat™. This catalogue was created using 
a large mosaic of observations with the VLT MUSE integral field spec- 
trograph and contains both LOS velocity measurements and metal- 
licities for more than 300,000 stars within the half-light radius of 
@ Cen. Although we could successfully cross-match five of the seven 
fast-moving stars, their signal-to-noise ratio is typically too low (about 2) 


for reliable velocity measurements, in particular for the four fastest, 
innermost stars. 

We could, however, obtain an LOS velocity value for star E (Vios = 
261.7 + 2.7 km s”) and star F (Vios = 232.5 + 4.0 km s”). These velocities 
are very close to the systemic LOS velocity of w Cen (232.99 + 0.06 kms”; 
ref. 13), confirming their membership in the cluster, as the Milky Way 
foreground is centred at Vios = O witha dispersion of 70 km s” (ref. 18). 
However, as the relative LOS velocity with respect to the cluster is low 
and we have only those two velocities for these outer stars, the LOS 
velocities do not add stronger constraints onthe IMBH. For this reason, 
we did not include them in the rest of our analysis. 


Testing the robustness of the assumed escape velocity 

Varying the parameters of the N-body models. Because we use 
the escape velocity of w Cen (assuming no IMBH is present) as the 
threshold for determining whether a star is considered fast or not, it 
is important to verify the robustness of the escape velocity value. We 
adopt an escape velocity Ves = 62 km s” (ref. 16); we have verified this 
value based on fitting similar N-body models to several state-of-the-art 
datasets, including MUSE LOS velocity dispersion measurements” 
and HST proper-motion-based dispersion measurements” for the 
central kinematics and Gaia DR3 measurements at larger radii using 
an assumed distance of 5.43 kpc. We varied both the assumed initial 
stellar mass function (using either the canonical Kroupa IMF™ or the 
bottom-light IMF derived in ref. 55) and the black hole retention frac- 
tion (assuming values of 10%, 30%, 50% or 100%). Despite changes to 
the central mass-to-light (M/L) ratio between the models, the central 
escape velocity changes only minimally, with a range of values from 
61.1 km s” to 64.8 km s”. Adopting any of these values leaves our sample 
of seven central stars above the escape velocity unchanged. 


An independent test using surface-brightness profiles. As a sec- 
ond test, independent of the N-body models, we calculated an escape 
velocity profile based on a surface brightness profile using various 
surface brightness profiles and dynamical models from the literature. 
We started by parameterizing the surface brightness profile using 
multi-Gaussian expansion (MGE)°** models. We then converted the 
surface brightness to a mass density using several mass-to-light ratios 
and distances used in the literature. From the mass density, we can 
derive the gravitational potential (@(r)). The escape velocity profile 
is then given by 


Vege (r) = [2(P (idav 0) - @(r,0)) 


with the tidal radius riqa = 48.6’ = 74.6 pc from refs. 57,58. 

These tests showed that the central escape velocity does not depend 
strongly on the stellar mass distribution in the centremost region, 
instead it is dominated by the global M/L ratio and the assumed dis- 
tance. The early dynamical models in the IMBH debate assumed both 
a distance of 4.8 kpc (ref. 59) and an M/L of 2.6 (vdMA10; ref. 7) or 2.7 
(NO8; refs. 3,4). With these values our tests give a central escape veloc- 
ity of 55.4 km s” (vdMA10) and 56.9 km s” (NO8). If we would also use 
the 4.8 kpc distance to scale the proper motions, this gives a cutoff of 
2.43 mas yr“ (vdMA10) and 2.48 mas yr” (NO8), close to the adopted 
cutoff at 2.41 mas yr’ and not changing the sample of the seven 
fast-moving stars detected. Owing to the parallax measurements of 
the Gaia satellite and updated kinematic distance measurements, the 
distance to w Cen was robustly redetermined and larger values have 
been found (5.24 + 0.11 kpc; ref. 60; 5.43 + 0.05; ref. 12). A dynamical 
model using the same surface brightness profile as ref. 3 but a larger 
distance of 5.14933 kpc was presented in ref. 8; this study found an 
M/L of 2.55*0 33. Using these values, the central escape velocity derived 
from the surface brightness profile is 61.1 km s” (equivalent to a proper 
motion of 2.51 mas yr”); again not changing our fast-moving stars sam- 
ple. Finally, varying the distance of any model by 0.2 kpc while holding 


the M/L constant results in variation in the escape velocity of the order 
of +3 kms”. These results show that our fast-moving star limit 
(Vese = 62 km sat a distance of 5.43 kpc) is consistent with the escape 
velocity values directly derived from surface-brightness profiles and 
several dynamically estimated M/L values. To visualize the escape veloc- 
ity, we calculated the escape velocity using the surface-brightness 
profile in ref. 3, an M/L of 2.4, and a distance of 5.43 kpc as found from 
the N-body models with a cluster of stellar mass black holes. The result- 
ing profile is shown in Extended Data Fig. 5f. The predicted escape 
velocity is flat out to approximately 50”, a property shared by all of the 
calculated escape velocity profiles. This makes the detection of the 
fast-moving stars only in the central few arcseconds more compelling. 


Anempirical confirmation of the central escape velocity. We make 
one final, and relatively model-independent, empirical confirmation of 
the central escape velocity based on the distribution of 2D velocities in 
the innermost region of w Cen (Extended Data Fig. 4). As we have seen 
in the analysis above, the escape velocity varies only slightly within 
the inner about 50” of the core of w Cen. Furthermore, the velocity 
dispersion profile is relatively flat in the innermost 10”, with a value 
of about 20 km s” (refs. 6,35,52). Therefore, we would expect rather 
similar distributions of stellar velocities in both the centre (0” < r < 3”) 
and an outer ring at (3” < r < 10”). Although we observe a clear excess 
of fast-moving stars in the inner 3 arcsec, there is a sharp cutoff very 
close to the adopted escape velocity in the (3” < r < 10”) bin. Although 
there is a total of 2,090 stars, there is only one star witha velocity sig- 
nificantly faster than the escape velocity (instead of 17 stars expected 
froma 2D Maxwell-Boltzmann distribution with dp = 20 km s”). This 
suggests that the stars with these velocities have escaped the central 
region. This outer fast-moving star has a 2D velocity of 75.8 kms ‘and 
is at a radius of r= 9.5”. From the density of Milky Way contaminants 
with apparent velocities above the escape velocity, we would expect 
about 0.7 foreground stars in the (3” < r < 10”) region; therefore, this 
fast-moving star is consistent with being a Milky Way foreground star. 


The escape velocity provides a minimum black hole mass 
The escape velocity for an isolated black hole is given by 


2GMpų 
Vesc,BH= \ oa 


In œ Cen, we have to take into account the potential of the globular 
cluster as well. If we assume this to be constant over the very small 
region in which we found the fast-moving stars (an assumption that 
agrees with the published surface brightness profiles, see Extended 
Data Fig. 6f), we obtain 


= |y2 2 
Vesc, total = döke + Vesc, cluster 


Ifa star at the distance of rsp with a velocity v,, is bound to the black 
hole, we can calculate the following lower limit on the black hole mass: 


2 2 2 2 
U3) Vo 
Men 7 I to) = Usciti 2, a) B Uku ciustéi 


A lower limit can also be calculated if the LOS velocity and distance 
are not known, aS Vsp > Vap and rp = Fop. AS we do not know the exact 
2D position of the black hole relative to the fast-moving stars, we cal- 
culated this lower limit for all stars and a grid of assumed 2D locations 
around the AvdM10 centreć. Each individual star alone would allow for 
a very low mass, as the location of the black hole could coincide with 
the star (Extended Data Fig. 5a-d). However, combining these limits 
for all stars gives a higher minimum black hole mass (Extended Data 
Fig. 5e). If we assume that all five robustly detected stars are bound to 
the black hole, the lower limit is about 8,200M, and the minimum mass 
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location is only 0.3” away from the AvdM10 (ref. 6) centre at the location 
RA = 201.6967370° and Dec = —47.4795066°. If we assume that the two 
most constraining stars are just random foreground contaminants, 
whichis ruled out at the 3a level (P = 0.0026), this limit drops to about 
4,100M,, still well within the IMBH range. 


Acceleration measurements 

The astrometric analysis in the catalogue” considered only linear 
motions of the stars. If there is a massive black hole present near the 
centre, we might also be able to measure accelerated motion of the 
closest stars, allowing for a direct mass measurement of the black hole. 
With an IMBH mass of 40,000M, and at a radius of 0.026 pc (1” on the 
sky), the acceleration of a star would be 0.25 km s” yr” (or 0.01 mas yr”). 
This is at the limit of the precision of our current dataset. With a 20-year 
baseline, we expect only a deviation of 0.05 pixel from a linear motion. 
For bright stars, the astrometric uncertainty can be as low as 0.01 pixel; 
however, for the faint fast-moving stars we have detected, the errors 
are significantly larger. 

To constrain these possible accelerations, we repeated the fit of the 
motion of each star enabling the addition of a quadratic component. 
The results for this fit are shown in Extended Data Table 1. The accelera- 
tions of the stars are consistent with zero within 30, but two stars have 
more than 2o acceleration measurements. The errors on our accelera- 
tion measurements lie between 0.004 mas yr’ and 0.03 mas yr? and 
are, therefore, of a magnitude similar to the expected acceleration 
signal. The strongest acceleration is shown by star B, which has been 
excluded from the robust subset of fast-moving stars because its proper 
motion is not 3g above the escape velocity. Owing to the proximity ofa 
bright neighbour star, we do not deem this acceleration measurement 
to be reliable. 

As the LOS distances to the fast-moving stars are unknown, it is not 
possible to place direct constraints on the IMBH mass using the upper 
limits on accelerations. If no acceleration is detected, as is the case 
for the centremost star A, this could mean that either the black hole 
is not very massive or the LOS distance of star A to the black hole is 
large. Combining the measurements for the ensemble of fast-moving 
stars and making some assumptions on their spatial distribution still 
enables us to use the acceleration limits to place further constraints 
on the black hole mass and its location. This is described in the next 
section. 


Markov chain Monte Carlo fitting of the acceleration data 
Assuming the fast-moving stars are bound to the IMBH, we can model 
the stars as being on Keplerian orbits around the IMBH. We used Bayes- 
iananalysis to sample the posterior distribution for the unknown mass 
and position of the black hole. In this analysis, there were eight free 
parameters: black hole mass; its on-sky x- and y-positions; and five 
LOS distances between the black hole and each fast-moving star. This 
analysis makes use of the available astrometric observations but stops 
short of modelling individual stellar orbits that would introduce more 
free parameters. 

We use a likelihood function with these eight free parameters and 
give the likelihood based on the observed on-sky x- and y-acceleration, 
proper motion and position of the five robustly measured fast-moving 
stars. For each star, we calculated a first likelihood term based on the 
modelled acceleration dpodeteq USING a Gaussian distribution with mean 
Aobservea aNd Width equal to the acceleration uncertainty. The second 
term in the likelihood accounts for the escape velocity constraints 
and is kept constant if the observed 2D velocity of the star is below 
the modelled escape velocity. For stars with 2D velocities above the 
modelled escape velocity, the likelihood is a Gaussian distribution 
with Mean yp — Vesctota aNd Width equal to the uncertainty in observed 
proper motion. 

We make these assumptions about the model: (1) The black hole mass 
is between 1M, and 100,000M.,, because a black hole mass beyond 


this upper limit is ruled out by our N-body models. (2) The black hole 
is located within the distribution of the fast-moving stars. We use a 
Gaussian prior in the black hole x- and y-positions with a mean equal 
to the mean position of the fast-moving stars and width equal to their 
one-dimensional positional standard deviation, 0,,,,, = 0.0221 pc; we 
also use a cutoff at +0.16 pc. (3) The stellar positions are isotropically 
distributed around the black hole. We model the LOS positions of the 
stars relative to the black hole using a Gaussian distribution with mean 
O and width o,,,,;- 

The posterior was sampled using a Markov chain Monte Carlo 
(MCMC) ensemble sampler implemented using the package emcee 
using recommended burn-in and autocorrelation corrections. We show 
the posterior distribution for the black hole mass in Extended Data 
Fig. 6a. The 99% confidence lower limit (21,100M.) is significantly 
higher than that derived from escape velocity constraints alone, 
whereas the upper limit on the mass is not well constrained. We also 
find a position for the black hole east of the AvdM10 centre, with 
Ax=-0.017'0-0" pc and 4y = 0.011709. pc (Extended Data Fig. 6b). 
The coordinates of the MCMC based centre estimate are RA = 
201.6970128° and Dec = —47.4794533°. We note that the black hole 
location estimate is dominated by the marginal 2g acceleration signal 
of star D, which is the faintest star in the sample; follow-up studies are 
required to obtain more precise acceleration measurements. 


N-body models 

Apart fromthe analysis of stars with velocities above the escape veloc- 
ity, we also used a set of existing N-body models with and without central 
IMBHs to get further constraints on the IMBH mass. We compared the 
simulations to the full velocity dispersion and surface density profile 
of w Cento determine the best-fitting model and the mass of acentral 
IMBH. The set of models and the details of the fitting procedure are 
described in detail in refs. 5,9. We note that these models have been 
presented already in the literature, but the fits to these models have 
been updated to incorporate the most recent Gaia DR3 data. 

Inshort, the models started from King profiles” with central concen- 
trations between c = 0.2 and c = 2.5 and initial half-mass radii between 
r,=2 pcandr, = 35 pc. Inthe models with an IMBH, we varied the mass 
of the IMBH so that it contains either 0.5%, 1%, 2% or 5% of the cluster 
mass at T= 12 Gyr when the simulations were stopped. The models 
with an IMBH assumed a retention fraction of stellar-mass black holes 
of 10%, whereas in the models without an IMBH we varied the assumed 
retention fraction of stellar-mass black holes between 10% and 100%. At 
the end of the simulations, we calculated surface density and velocity 
dispersion profiles for each N-body model and then determined the 
best-fitting model by interpolation in our grid of models and using x? 
minimization against the observed velocity and surface density profile 
ofw Cen. 

The velocity distributions from observations and the models are 
shown in Extended Data Fig. 7. We compare the distribution of meas- 
ured 2D stellar velocities in the inner 10” of w Cen with the various 
models using a Kolmogorov-Smirnov test. Moreover, we compare 
the fraction of fast-moving stars in the innermost 3 arcsec (Extended 
Data Table 2). Models without an IMBH and witha 20,000M, IMBH are 
both strongly excluded by both the overall velocity distribution and the 
complete lack of fast-moving stars. The overall velocity distribution is 
in best agreement with the 47,000M, distribution, whereas the frac- 
tion of fast-moving stars is best matched by the 39,000M, simulation. 
This tension might be alleviated in future models that contain both an 
IMBH and a cluster of stellar mass black holes. We caution that these 
simulations have smaller numbers of stars than observed and that 
there can be substantial variations in the distribution of central stars 
due to strong encounters with remnants and binaries. Nonetheless, 
the simulations suggest that black holes with masses of M < 50,000M, 
are consistent with the observed distribution of central velocities and 
fast-moving stars, whereas the no-IMBH case and significantly more 


massive black holes are disfavoured because of an overprediction 
of fast-moving stars. Updated N-body models fit to the OoMEGACat 
kinematic data and dynamical modelling of these same datasets with 
Jeans models are currently underway. 


Comparison with S-stars in the Galactic centre 

The black hole indicated by our fast-moving star detection is only the 
second after Sgr A*, for which we can study the motion of multiple 
individual bound stellar companions. Therefore, the extensively meas- 
ured stars around Sgr A* provide a unique comparison point to our 
fast-moving star sample. We compare the motions of the stars in the 
S-star catalogue in ref. 17 with our w Cen fast-moving star sample in 
Extended Data Fig. 8. When taking into account the different distances 
and the approximate black hole mass ratio of 100, the motions show 
similar amplitudes. However, the density of tracers in w Cen is signifi- 
cantly lower, despite the greater depth of the observations. 


Data availability 


The data used in this paper are based on archival observations taken 
with the HST that are freely available in the Mikulski Archive for 
Space Telescopes. All used observations have been grouped under 
a DOI (https://doi.org/10.17909/26QJ-GO90)*. Moreover, the full 
proper-motion catalogue is made public along with the respective 
publication”. 


Code availability 


We used the following Python packages to perform the analysis: mat- 
plotlib™, scipy®, numpy®, astropy” and emcee“. The N-body simula- 
tions were run with the publicly available NBODY6 code®’. We can share 
the code used in the data analysis upon request. 
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Extended Data Fig. 1| Number density of fast-moving stars. The black rise of the density towards the AvdM10° centre. Based on the number density 
markers with errorbars (10; based on Poisson statistics) show the measured of fast-moving stars at large radii (whichis consistent with predictions froma 
number density of robustly detected fast-moving stars determined in radial Besancon Milky Way model, dashed line), only ~0.073 fast-moving stars are 
bins with respect to w Cen’s center. A constant density of Galactic fore-/ expected within the central 3” compared tothe 5 observed stars. 
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Extended Data Fig. 2| Astrometry and multi epoch imaging for all fast- 
moving stars. Each row (a-g) shows the astrometry and imaging for one 

of the seven fast-moving stars. The left column shows the raw astrometric 
measurements used to determine the proper motions, colour-coded by the 
epoch of their observation. The centre column shows the linear and quadratic 
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fits to both the R.A. and Dec. position change of the stars. Errorbars correspond 
to the loerror on the individual position measurements. The right column 
shows stacked images from 2002 (ACS/WFC F625W) and 2023 (WFC3/UVIS 
F606W); the positions of the fast stars are marked witha pink open circle. 
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Extended Data Fig. 3 | Completeness of the catalogue and photometric 
diagnostics. Panel a compares the completeness of the various available proper 
motion datasets for the core of w Cen, using histograms of the magnitude 
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here. The grey dots in b-d show the mean of photometric diagnostics for the 
raw PSF photometry measurements for the fast-moving stars compared with 
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the bulk of stars in the catalogue. The first panel shows the QFIT parameter, 
given by the linear correlation function between the PSF and the pixel values 
inthe image. The second panel shows the radial excess parameter, a parameter 
that compares the residual flux inside versus outside of the fit aperture. 

The third panel shows the flux ratio between the flux of a star itself and its 
neighbouring sources. All five robustly measured stars show typical behaviour 
for their magnitude, while the excluded fast-moving stars (B, G) areinfluenced 
by bright neighbouring sources. 
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Extended Data Fig. 4 | Empirical verification of the escape velocity. Panels a 
and b show histograms of the observed 2D velocity distribution in the very 
centre (O” <r <3”; c) andinan outer ring (3” <r < 10”; d). While the lower velocities 
are well described by a 2D Maxwell-Boltzmann distribution with op =20 km s” 
(marked witha solid black line, the dashed black lines refer to alternative 
distributions with o,)=17 kms ‘and dp =23 kms’), there are clearly notable 
differences at higher velocities. Those become especially visible in the 
cumulative normalized histogram shown in d and the zoom-in ine: while the 
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distribution between (0” <r < 3”, blue line) shows an excess of fast-moving stars, 
the distribution at larger radii (3” <r < 10”, orange line) shows a clear deficit of 
stars at velocities larger than the escape velocity, making the used escape 
velocity threshold very plausible. Even though the sample is ten times larger, 
there is only a single star witha velocity significantly larger than Vese- This star 
has a2D velocity of 75.8 kms ‘and is at a radius of r = 9.5”. Itis consistent with 
being a Milky Way foreground star. 
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escape velocity of the cluster indicates that they are bound to a massive object. based escape velocity profile in blue, using the surface brightness profile from’ 
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(a-d) show the contours of the minimum black hole mass indicated by the four The radii are measured with respect to the minimum mass centre shown ine. 
centremost robustly measured fast-moving stars. By combining the minimum The shaded regions indicate the uncertainty introduced by an assumed error of 
black hole mass constraints from each of the fast-moving stars, we can find the +0.2 kpc onthe distance. The surface-brightness profile based escape velocity 
position that allows for the lowest IMBH mass (e). This analysis indicates a firm is compatible with the adopted value of Ves- = 62 kms. The profile without an 
lower limit of around 8,200 M.. Our minimum mass location only differs by IMBHis also nearly flat in the inner ~50”, justifying the assumption of a flat 
~0.3 arcsec from the AvdM10° centre. This result does not significantly change profile in the innermost region. 
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contours shown correspond to the 1-, 2-, and 3-sigma levels of the distribution. distribution from an MCMC run including stars B and G. 
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Extended Data Fig. 7| Comparison of the observed velocity distribution 
with N-Body models. a, 2D velocity distribution for the stars in the inner 
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without an IMBH predicts no stars above the escape velocity; the 47,000 M, 
model predicts anumber close to our observations. b, Comparison of the 
normalized, cumulative distribution of stellar velocities for our data and five 
different N-body models. 
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Extended Data Fig. 8 |Comparison with the Galactic Centre. In this figure 
we compare the observed physical motion of our fast star sample with the 
stars orbiting the black hole Sgr A* in the Galactic Center”. The physical scale 
probed by the fast-moving stars is similar to that probed by the S stars inthe 
Milky Way centre; however, the density of these tracers is lower. Owing to the 
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Extended Data Table 1 | Detailed astrometric and photometric information of the fast-moving stars 


a 
General Properties 
Star Catalog ID R.A. Dec. Fproj. Nused 
[degree] [degree] [J 
A 532867 201.6967263 -47.4795835 0.265 308 
B 475236 201.6968888 -47.4797138 0.537 353 
C 1422379 201.6966378 -47.4793672 0.870 262 
D 476492 201.6968346 -47.4793233 0.886 195 
E 509644 201.6973080 -47.4797545 1.333 430 
F 476467 201.6977125 -47.4792625 2.408 427 
G 510061 201.6961340 -47.4790585 2.506 338 
b 
Velocity and Acceleration Measurements 
Star PM R.A. pacos 5 PM Dec. ps Total PM V2D Acceleration R.A. Acceleration Dec. 
mas yr? mas yr? mas yri km s^ mas yr- mas yr- 
A 3.563+0.038 2.564+0.055 4.390+0.044 113.0+1.1 -0.0069+0.0083 (0.80) 0.0085+0.0098 (0.90) 
B 2.167+0.182 1.415+0.081 2.588+0.159 66.6+4.1 0.0702+0.0239 (2.90) 0.0228+0.0157 (1.50) 
C 1.117+0.127 3.514+0.056 3.687+0.066 94.9+1.7 0.0028+0.0333 (0.10) -0.0060+0.0123 (0.50) 
D 2.559+0.082 -1.617+0.061 3.027+0.076 77.9+2.0 0.0357+0.0177 (2.00) -0.0194+0.0162 (1.20) 
E -2.149+0.025 1.638+0.037 2.702+0.030 69.6+0.8 0.0072+0.0042 (1.70) -0.0009+0.0075 (0.10) 
F 0.436+0.017 -2.584+0.016 2.620+0.016 67.4+0.4 0.0052+0.0038 (1.40) -0.0015+0.0038 (0.40) 
G -1.317+0.098 2.207+0.062 2.571+0.073 66.2+1.9 -0.0197+0.0267 (0.70) 0.0173+0.0170 (1.00) 
c 
Photometric Properties 
Star mrcosw QFITrFeoew RADX Feéosw f/f. Feoew MF814W QFITF814w RADX Fs14w f/f. F814w 
A 23.373+0.009 0.976 (p =0.716) -0.006 (p =0.192) 0.366 (p =0.532)  22.134+0.027 0.945 (p =0.170) -0.004 (p =0.249) 0.363 (p =0.537) 
B 22.778+0.014 0.940 (p =0.096) 0.082 (p =0.981) 0.122 (p =0.322) 21.625+0.026 0.956 (p =0.128) 0.072 (p =0.964) 0.494 (p =0.639) 
C 23.630+0.009 0.940 (p =0.283) 0.050 (p=0.8249 0.147 (p=0.320) 22.342+0.034 0.916 (p=0.109) -0.014 (p =0.136) 0.327 (p =0.514) 
D 24.10840.017 0.952 (p =0.585) -0.008 (p =0.216) 0.064 (p =0.180) 22.745+0.037 0.952 (p=0.291) -0.056 (p =0.022) 0.079 (p =0.198) 
E 20.112+0.004 0.996 (p =0.131) 0.014 (p =0.950) 0.424 (p =0.923) 19.373+0.009 0.995 (p =0.109) 0.018 (p =0.938) 0.622 (p =0.975) 
F 21.348+0.004 0.997 (p =0.635) -0.002 (p =0.340) 0.028 (p =0.253) 20.391+0.004 0.998 (p =0.799) 0.004 (p =0.587) 0.032 (p =0.232) 
G 20.888+0.009 0.980 (p =0.031 -0.003 (p =0.219) 2.049 (p =0.997) 20.069+0.014 0.983 (p =0.053) 0.007 (p =0.730 0.971 (p =0.984) 


The table in a shows the position, the number of astrometric measurements N,,.eg, and the projected distance from the AvdM10 centre for each of the fast-moving stars. The table in b lists the 
individual proper motion components, the total proper motion, the inferred 2D velocity and the measurements of the acceleration for the seven fast-moving stars. All shown errors correspond 
to the 10 errors, which were estimated by scaling the formal errors on the parameters by dxa. of the respective fit. The strongest accelerations are shown by Star B; however, this star has been 
discarded from the set of robustly measured stars because of large astrometric errors. All robust stars show an acceleration consistent with zero to within 20. Finally, the table in ¢ lists several 
photometric properties of the fast-moving stars in two filters, including the measured brightness and several photometric diagnostics (the quality of fit (QFIT) parameter, the radial excess 
(RADX) parameter and the flux ratio between the flux of each source and the flux of neighbouring sources). Together with the photometric diagnostic values we show their percentile with 
respect to stars with a similar magnitude (p). 


Extended Data Table 2 | Comparison of the observed velocity distribution with N-body models 


IMBH mass 

fraction in Mimsn KS test N3” fast 
model Mo Nio”tota p-Value N3” fast  N3rtota  N3”total 
Observed Data 2,324 - 7 218 3.2% 
(0) 0 1,148 6.21e-7 0 113 0% 
0.5% 20,000 1,342 8.60e-10 0 113 0% 
1.0% 39,000 868 0.0027 3 89 3.4% 
1.2% 47,000 901 0.27 7 85 8.3% 
2.0% 78,800 843 0.18 7 70 10.0% 


In this table we report results of the comparison of the observed distribution of 2D velocities with different N-body models (see also Extended Data Fig. 8). The first two columns indicate the 
relative and absolute mass of the IMBH in the N-body models. The third column shows the total number of stars within 10” of the cluster centre. In the fourth column, we show the results of 
a Kolmogorov-Smirnov test comparing the measured velocity distribution with the different N-body models. Finally, the last three columns compare the absolute number and the fraction of 
fast-moving stars in the centremost 3 arcseconds. 


